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^ Abstract 

P^P The recent LHC discovery of a Higgs-like resonance at 126 GeV suggests that the minimal 

supersymmetric standard model must be modified in order to preserve naturalness. A simple 
extension is to include a singlet superfield and consider the scale-invariant NMSSM, whose 
renormalizable superpotential contains no dimensionful parameters. This extension not only 
solves the /i-problem, but can easily accommodate a 126 GeV Higgs. We study the naturalness 
Qh of the scale-invariant NMSSM taking into account the recent constraints from LHC searches, 

flavor physics and electroweak precision tests. We show that TeV-scale stop masses are still 
allowed in much of the parameter space with 5% tuning for a low messenger scale of 20 TeV, 

i-C split families (with third-generation sleptons decoupled) and Higgs-singlet coupling A of order 

one. For larger values of the Higgs-singlet coupling, which can relieve the tuning in the Higgs 

fN| VEV, an additional tuning in the Higgs mass limits increasing the (lightest) stop mass beyond 

1.2 TeV, the gluino mass above 3 TeV, and electroweak charginos and neutralinos beyond 400 
GeV for a combined tuning better than 5%. This implies that the natural region of parameter 

^ space for the scale-invariant NMSSM will be fully explored at the 14 TeV LHC. 
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1 Introduction 



The recent discovery of a resonance with Higgs-like properties and mass around 126 GeV at the 
Large Hadron Collider (LHC) provides compelling evidence that the Higgs mechanism is respon- 
sible for the generation of elementary particle masses in the Standard Model (SM). An important 
next question is to ascertain the nature of this new resonance by measuring its properties and cou- 
plings. Preliminary indications already suggest an enhancement of the Higgs decay rate to photons 
which would signal the presence of new light states beyond the standard model. A well-motivated 
possibility is to identify the new resonance with the lightest Higgs boson in the supersymmetric 
standard model. Supersymmetry (SUSY) provides a natural solution to the hierarchy problem while 
simultaneously allowing for gauge coupling unification and a suitable dark matter candidate [ ]. 
However, a Higgs boson at 126 GeV is in extreme tension with maintaining naturalness in the 
Minimal Supersymmetric Standard Model (MSSM). As is well-known, large radiative corrections 
from stop loops are required to enhance the Higgs quartic coupling which originates at tree-level 
from D-terms. This is accomplished by either having large (S> 1 TeV) stop masses, in the case 
of negligible mixing in the stop sector, or a large trilinear Af-term. In either case, a fine-tuned 
cancellation caused by this little hierarchy is needed to obtain the correct electroweak symmetry 
breaking scale. 

To alleviate the dependence on large stop masses or the ^4^-term in obtaining a sufficiently large 
Higgs quartic coupling, the MSSM can be extended to include a gauge-singlet chiral superfield 
S in addition to the MSSM particle content. In this Next-to-Minimal Supersymmetric Standard 
Model (NMSSM; see [2] for a review) the singlet superfield couples to the up and down Higgs 
superfields, H u and H^, via the superpotential term XSH u Hd, where A is a dimensionless coupling. 
This modifies the upper bound on the Higgs mass at tree-level with respect to the MSSM such that 

m 2 h < m| cos 2 2(3 + AV sin 2 2/3, (1.1) 

where v is the electroweak vacuum expectation value (VEV) and tan f3 denotes the ratio of the up 
and down Higgs VEVs. The extra tree- level contribution in (1.1) from the singlet coupling A is 
maximal for tan /3 ~ 1 and can thus increase the Higgs mass, thereby alleviating the need for large 
stop masses or a large ^-term and the resulting fine-tuning problem. In addition, values of A > 1 
further ameliorate the fine-tuning problem by suppressing the sensitivity of the Higgs VEV with 
respect to changes in the soft masses [3; 4]. 

In addition to the Higgs-singlet coupling there can also be other singlet field couplings, including 
explicit mass terms in the superpotential. These mass terms can be avoided by considering the 
scale-invariant version of the NMSSM which is characterized by forbidding all dimensionful terms in 
the superpotential at the renormalizable level. This can be accomplished via a discrete Z3 symmetry 
under which all superfields are charged, only allowing the superpotential terms \SH u Hd + (k/3)S 3 
(the resulting domain wall problem in the early universe can be addressed as discussed e.g. in [5; 6]). 
This scenario has the additional virtue of providing a solution to the /^-problem encountered in the 
MSSM once the scalar component of S acquires a VEV of order the electroweak scale. 

In this paper, we study the parameter space of the scale-invariant NMSSM. We search for 
regions that can incorporate a Higgs at 126 GeV and that are in accordance with all experimental 
constraints. Our main focus are the constraints on the sparticle spectrum coming from naturalness. 
To quantify the naturalness, we use the usual logarithmic-derivative measure of fine-tuning. We 
not only include the one-loop radiative corrections coming from supersymmetric running in the 
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calculation of the fine-tuning measure but also contributions from the Coleman- Weinberg potential, 
which are particularly important for (s)particles coupling strongly to the Higgs-singlet sector. In our 
analysis we will make a number of assumptions. First, the only superpartners that we consider in 
addition to those in the Higgs-singlet sector are stops, sbottoms and gauginos. All other sfermions 
are taken heavy enough to be decoupled from our effective theory. Decoupling these sfermions, 
corresponding to split families (including the third generation sleptons), helps to satisfy several 
collider and flavor constraints but does not affect the naturalness of the Higgs sector because these 
sparticles couple only weakly to the Higgs. Furthermore, in order to allow for a large value of A 
and to minimize loop corrections to the soft parameters, we consider a low cutoff of our effective 
theory, A mess ~ 20-1000 TeV. This scale could be associated with a messenger sector like in gauge 
mediation. Alternatively, it could correspond to a cutoff below which the Higgs sector emerges 
as composites of an underlying strongly coupled theory as in accidental SUSY [7-9] or fat Higgs 
models [10]. Finally, we allow all values for A and k that are consistent with perturbativity up to 
the cutoff scale A mess - This covers not only the parameter space of the usual NMSSM (for which 
one typically requires perturbativity of the couplings up to the GUT scale) but also includes values 
A > 0.7 that develop a Landau pole below the GUT scale as considered in ASUSY models [3]. As 
shown recently [11], a Landau pole in the coupling A is not necessarily detrimental to maintaining 
gauge coupling unification. 

The assumptions that we make - split families (with third-generation sleptons decoupled), a 
low cutoff A mcss and a possibly large coupling A - correspond to an optimal-case scenario for the 
(scale- invariant) NMSSM from the viewpoint of naturalness. Indeed, if the families are not split, 
stops and sbottoms would tend to be much heavier due to the strong collider constraints on first- 
generation squarks. Similarly, a larger cutoff scale A mess would result in larger loop corrections to 
soft parameters from supersymmetric running. Large values A > 0.7, on the other hand, alleviate 
the fine-tuning problem as discussed below Eq. (1.1). This ensures that our constraints from 
naturalness are conservative in the sense that a different implementation of the (scale-invariant) 
NMSSM would have a smaller region of parameter space for a given amount of fine-tuning. 

The amelioration of the Higgs VEV tuning at large A seems to allow for stop masses above the 
TeV scale for the same level of tuning as in the MSSM. However, we point out that there is also 
a competing effect. The bound in Eq. (1.1) implies that the Higgs mass at tree- level grows with A 
if tan /3 is small. Electroweak precision tests, on the other hand, require that tan /3 ~ 1 for large 
A. This means that the Higgs mass typically becomes larger than 126 GeV in the region A > 1. A 
tuning in the various contributions to the Higgs mass is then required to bring the Higgs mass down 
to 126 GeV. This additional tuning prevents the stop masses to be raised much above the TeV 
scale if the theory is supposed to remain natural. The two competing effects can be incorporated 
into a single measure of tuning, which we define as the product of the Higgs VEV and Higgs mass 
tunings. We find that this combined tuning measure is minimized for A ~ 1, leading to restrictions 
on a natural sparticle spectrum. 

We use the numerical package NMSSMTools 3.2.1 to compute the sparticle spectrum, but mod- 
ify the program to allow for a small messenger scale A mess and a large coupling A. We impose 
constraints from direct searches at colliders and flavor physics. Several such constraints are already 
implemented in NMSSMTools. We complement these with the recent limits on heavy CP-even 
Higgses from ATLAS and CMS, as well as limits from stop and gluino searches. In addition, we 
require that the Peskin-Takeuchi parameters S and T are consistent with their measured ranges. 
Finally, we verify that the LSP relic abundance does not overdose the universe and impose con- 
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straints from direct detection. We then perform a numerical scan of the parameter space of the 
scale-invariant NMSSM which satisfies these constraints and optimize the scan for regions that have 
a small tuning in the Higgs VEV and a Higgs mass near 126 GeV. The results from the scan leads 
to restrictions on a natural sparticle spectrum. It is found that for a total tuning which is better 
than 5%, and A mess ~ 20 TeV, the lightest stop must have a mass below 1.2 TeV, the gluino should 
be below 3 TeV, and the electroweak charginos and neutralinos must lie below 400 GeV. These 
mass ranges still leave room for naturalness to be explored by the 14 TeV LHC. Furthermore if the 
total tuning is relaxed to be better than 1% then are a few points where a large enhancement in 
the Higgs diphoton signal strength can be obtained. 

There have been a number of previous analyses [3; 4; 12] where non-minimal supersymmetry 
with large A was studied in the context of naturalness. These differ from our work in that they 
do not perform a full scan of the parameter space. Furthermore, Refs. [3; 4] include explicit mass 
terms (such as a /x-term and a singlet supersymmetric mass) in the superpotential. These papers 
also neglect the quantum effects accounted for by the Coleman- Weinberg potential which we find 
to be important. Naturalness and a 126 GeV Higgs boson in the context of the NMSSM with 
perturbativity up to the GUT scale (i.e. A < 0.7) was considered e.g. in [13-16]. Our work also 
complements the recent analysis in Ref. [17], where stop masses were restricted to lie below 500 
GeV. 

The paper is organized as follows. In sec. 2, we present the Higgs sector of the scale-invariant 
NMSSM and include the contributions from the effective potential. A brief description of the 
numerical scan with the range of the model parameters is given in sec. 3. The naturalness constraints 
are discussed in sec. 4, where fine-tuning can arise from both the electroweak VEV and the Higgs 
mass. Requiring small combined tuning leads to constraints on the stop masses and the Higgs- 
singlet coupling A. Further constraints are described in sec. 5, which arise from collider limits, 
electroweak precision tests, flavor physics and dark matter. In addition in sec. 6 we also present 
predictions relevant for Higgs searches and comment briefly on the consequences for SUSY searches. 
Our conclusion is given in sec. 7 and in the Appendices we present, for completeness, relevant 
expressions for sparticle mass matrices and the effective potential used in our analysis. 

2 The scale-invariant NMSSM 

We consider a simple extension of the MSSM, where an extra chiral singlet superfield S is added 
to the usual MSSM field content. The superpotential couplings of the singlet are constrained by 
a Z3 symmetry which forbids linear and quadratic terms. This leads to the superpotential of the 
scale-invariant NMSSM ^ 

Wnmssm = ^SH u Hd + -S 3 , (2.1) 

where A and k are dimensionless couplings. Including contributions from soft SUSY breaking and 
the D-terms, the resulting scalar potential reads 

V = (m 2 Hd + X 2 S 2 )\H d \ 2 + (m 2 Hu + X 2 S 2 )\H U \ 2 + X 2 \H d H u \ 2 + m 2 \S\ 2 + k 2 \S\ a 

+ [(a x S + \KS 2 )H u H d + ^S 3 + h.c] + V D , (2.2) 
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where mH u ^H d ^s^ a \ and a K are soft mass parameters. Furthermore, Vd denotes the D-term 
contribution to the potential, 

Vd = V U(1)Y + V SU(2)W = ^±M2 (|ff u | 2 - |P d | 2 ) 2 + | \H\H U \\ (2.3) 

where g± is the U(l)y and 52 is the SU(2) gauge coupling. To ensure CP-conserving VEVs, we shall 
assume that the Lagrangian of the Higgs sector is CP-invariant. By field redefinitions, it is then 
possible to choose X,a K ,v u and Vd real and positive, while the rest of the Higgs-sector parameters 
are real but can have both signs. For completeness, let us mention that the Yukawa couplings 
yt > and y\, > can similarly be chosen real and positive. 

The soft SUSY breaking terms are generated at a scale A mess , below which they receive loop 
corrections via renormalization group (RG) running. The scale A mess could be associated with the 
mass of messenger fields which communicate SUSY breaking from a hidden sector to the visible 
sector. Alternatively, A mess could be a cutoff below which NMSSM fields emerge from an underlying 
strongly coupled theory as, for example in the models of [7-10]. For definiteness, we shall refer to 
A mess as the messenger scale in the following. Our aim is to identify a region of parameter space 
with minimal fine-tuning for the scale-invariant NMSSM. In order to limit loop corrections to soft 
masses from RG running, we shall therefore focus on a relatively low messenger scale A mess = 20 TeV 
(though we also comment on the cases A mcss = 100 TeV and 1000 TeV). 

The NMSSM provides an additional tree- level contribution to the Higgs mass which is enhanced 
for large values of A (see below). Due to the fast running of A, values A > 0.7 at the electroweak 
scale lead to a Landau pole below the GUT scale. In order to accommodate larger A and thus a 
sizable tree-level contribution to the Higgs mass, we will only require perturbativity of the couplings 
up to the messenger scale. This is in the spirit of models of ASUSY [3; 4] and allows for values 
A < 2.3 at the electroweak scale. Note that some form of UV completion is required to kick in 
before the Landau pole in A (or any of the other couplings). The cutoff associated with this UV 
completion would have to coincide with the messenger scale A mess if a Landau pole occurs shortly 
above A mcss (e.g. for A « 2.3) but otherwise this is not necessarily the case. 

In order to account for quantum fluctuations we make use of effective field theory in the form 
of effective action methods that are based on the run-and-match procedure by which all relevant 
parameters are run with the RG scale, for scales greater than the masses involved in the theory. 
When the RG scale goes through a particular mass threshold, heavy fields decouple, leaving thresh- 
old effects that are used to match the effective field theory above and below the mass threshold. 
The effective action for the Higgs scalars, through a loop expansion, can be written as: 



„ ( 00 00 ^ 

s eff = y E w!^ - E Vn \ > ( 2 - 4 ) 



v n=0 n=0 

where (pi = H u ,Hd,S, the Zf account for wave-function renormalization, and the loop calculation 
involves a mass parameter \i r known as the renormalization scale. The well-known Coleman- 
Weinberg [18] expression for the one- loop effective potential in the DR scheme is given by 



M 2 \ 3 



log, ^y 2 



(2.5) 



where the trace runs over all particles in the effective theory, M denotes their masses and all 
dimensionless and dimensionful couplings of the theory are evaluated at [i r . The dominant contri- 
butions to V\ can be found in the Appendix. A sensible choice for the renormalization scale fi r that 
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makes the approximation more scale-independent and minimizes large logarithms arising from the 
top/stop sector is fi r = m so f t = ^/niQ 3 m U3 , where mg 3 and m U3 are the soft-breaking masses of the 
stop sector. Therefore, above the scale m so ft we have an effective supersymmetric theory where all 
couplings can be determined at any scale m so ft < Q < A mess using supersymmetric RG running. 
In this way, potentially divergent logarithms that involve the messenger scale A mess are taken into 
account by integrating the SUSY RG equations, while in principle all finite terms are accounted 
for by the effective action. Given the potentially large values of the couplings yt, A, k at the SUSY 
scale, we expect that the major contributions to the effective action come from the stop/top and 
Higgs-singlet sector. 

Minimizing the effective potential (including I-loop corrections) with respect to the neutral 
components of H u , and S, we find that the following minimization conditions must be satisfied 

m 2 Hu + X 2 v 2 s -(axvs + XKV 2 s )^^- r ^cos2^ + X 2 v 2 cos 2 ^ = 0, (2.6) 

9 

™L + ^ v s - («a^5 + XkvI) tan p H f- cos 2(3 + X 2 v 2 sin 2 (3 = 0, (2.7) 

d 2 

2 

m 2 s -XKv 2 sin2p + 2K 2 v 2 s + X 2 v 2 - s\n2(3 + a K v s = 0, (2.8) 

2vs 



where 



™H n = rn Hu + m Ha = m Hd + -^V u fh s = m s + -^V h (2.9) 

w d S 

and v 2 = v 2 + v^, tan/3 = v u /vd, m 2 z = \{g\ + g^)v 2 and all parameters are evaluated at the 
renormalization scale fj, r . From the first two equations, we can derive the following relations: 



Xv = 2 ■ ofl ~ m Hu - m H , - 2A v s , (2.10) 



aJ n (ax VS + ^V 2 s ) ^ ~2 , )N L>.,2 

sin 2/8 Hu 

^2 ^2 

z ~~ — cos 2/3 Hu Hd s ' ' ' 

Notice that the second condition (2.11) is similar to the one obtained in the MSSM if we iden- 
tify [i = Xvs and neglect Higgs-singlet contributions to the effective action. As can be seen from 
Eq. (2.4), the fields H u , H d and S do not yet have properly normalized kinetic terms. Absorbing 
the wavefunction renormalization factors Zjj u , Zn d and Z$ into the fields, the VEVs of the prop- 
erly normalized fields are related to the VEVs v u , Vd and vs found from (2.6)-(2.8) according to 
Wu,norm. = v u yj Zjj u and similarly for Vd and vs. In the following, we will drop the subscript 'norm.' 
and v u , Vd and vs will instead refer to the VEVs of the properly normalized fields. We then require 
that v = 174 GeV. 

Expanding the neutral components of the Higgs fields around their VEVs, 

o h u + ih U}I o h d + ih d ,i s + is/ 

" v 7 2 = ' H d = Vd+ ~7f ' 5 y/2 ' 

where the subscript / denotes the imaginary part, one finds the mass matrix of CP-even states: 

(M 2 \ - ^CP-even) ij + ^ 
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Here M^j p _ even denotes the tree-level contribution which in the basis (h u ,hd,s) is given by the 
matrix 

(m% sin 2 /3 + m 2 A cos 2 /3 (A 2 w 2 — | (m% + m\)) sin 2/3 2A/iu sin /3 — a\v cos /? — 2k/i« cos/? \ 
sin 2 /3 + m 2 z cos 2 /? 2A/iu cos /3 — a^u sin /3 ~ 2nfiv sin/3 J , (2.14) 

A 2 « 2 +m| + 6f^K 2 -AKW 2 sin2/? + 2a fe ^ / 

We adopt the convention that dotted entries are obtained from the fact that the matrix is symmetric 
and, for convenience, we have defined fi = \v$ and 

2 

sin2(3^ axVS + XKV ^ ' ^ 2 ' 15 ^ 

Furthermore, AM 2 denotes contributions from loop corrections to the effective potential which are 
given by 



The eigenvalues m 2 ([i r ) of the mass matrix (2.13) are running masses evaluated at the renor- 
malization scale fj, r . From these, the physical pole masses are obtained by solving the pole equation 

det [p 2 <% - (mf(fir) + AU t (p 2 ))] = , (2.17) 

where AII^p 2 ) = U^p 2 ) - 1^(0) and Ui(p 2 ) are the L>i?-renormalized self energies. 

It is useful to consider a basis in which the electroweak VEV is associated with only one linear 
combination of states hi '. To this end, we rotate 



sin /3 


cos /3 






( K \ 


— cos p 


sin/3 
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1 } 







(2.18) 

In this basis, the diagonal element of the tree- level mass matrix associated with h' becomes 

m h>h> = m z cos2 W + A 2 w 2 sin 2 2/3 . (2.19) 

This constitutes an upper bound on the lightest eigenvalue of the CP-even mass matrix at tree- 
level. The second term, which is not present in the MSSM, implies that the Higgs mass can be 
larger than the Z-boson mass already at tree level. In comparison to the MSSM, this alleviates the 
need for heavy stops to raise the Higgs mass to 126 GeV. 

Note that the bound is saturated (so that Eq. (2.19) is the Higgs mass) if the lightest CP-even 
state is purely h! . Since in this basis, the VEV is only carried by h' , such a Higgs would couple 
to SM fields just like the SM Higgs. Given the uncertainties in the production cross sections and 
branching fractions of the Higgs (and the hint for an excess in h — > 77), a non- vanishing admixture 
of h' with the orthogonal state H' and the singlet s is still permissible and typically present. This 
reduces 1 the Higgs mass relative to Eq. (2.19). In addition, loop corrections may either increase or 
reduce the Higgs mass. 



1 The state identified with the Higgs might not be the lightest CP-even state (see, for example [17]). In this case, 
the admixture raises the Higgs mass. We do not find any such points in our numerical scan and will therefore not 
consider this case further. 
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tan/3 


tan/3 > 0.08 


mQ 3 


Ag7Ti,Q 3 < mg 3 < 5 TeV 


Mi 


< Mi < 8 TeV 


A* 


< ITeV 


m U3 


Agin U3 < m U3 < 5 TeV 


M 2 


< M 2 < 8 TeV 


A 


< A < 3 




< m d . 3 < 8 TeV 


M 3 


0.5 TeV < M 3 < 8 TeV 




|k| < 2.75 


A t 


\A s A t \ < \A t \ < 5 TeV 


V 


174 GeV 




|A A | < 2TeV 


A b 


\A b \ < 8 TeV 


■A-mess 


20, 100, 1000 TeV 


A K 


\A K \ < ITeV 











Table 1: The ranges of the input parameters (defined at the renormalization scale m so ft) which we 
use in the numerical scan. Here A\ = a\/X, A K = a K /n and Ag£ denotes the one-loop contribution 
of gluinos to the soft mass parameters £ = At,mQ 3 ,m U3 . We choose a K > without loss of 
generality as discussed in sec. 2. 

3 Numerical Analysis 

Before proceeding to the discussion of the imposed bounds and results in detail, we give a brief 
description of how the numerical scan is performed. Readers who are mainly interested in our 
results could skip this section. 

As motivated in the Introduction, we decouple all sparticles from the first and second generation 
as well as third-generation sleptons. The remaining parameter space of the scale-invariant NMSSM 
is characterized by the 16 parameters given in Tab. 1. Note that we use the VEVs v u , and vg 
(or v, tan/3 and fj,) as input parameters instead of the soft masses w 2 ^, m # d an d rn 2 s . The latter 
can be obtained from the former by means of the minimization conditions (2.6) - (2.8). After fixing 
the Higgs VEV at v = 174 GeV, the messenger scale at A mess = 20 TeV (and all SM masses and 
couplings at their measured values), we perform a numerical analysis by sampling points in the 
remaining 14-dimensional parameter space defined at the renormalization scale m so ft = ^JmQ 3 m U3 . 
To this end, we employ a Markov Chain Monte Carlo using the Metropolis-Hastings algorithm with 
simulated annealing. We scan linearly in all parameters, i.e. we prefer large values in the scan, 
which helps to find the boundary of the region with small fine-tuning. The likelihood function in 
the scan is given by the product of a Gaussian for the Higgs mass (i.e. the state that we identify 
with the recently discovered resonance) centered at 126 GeV and a Gaussian for the fine-tuning 
centered at 0, i.e. the sampler prefers regions which have a Higgs mass close to the experimentally 
observed value and small fine-tuning in the electroweak VEV. 2 We will discuss the exact definition 
of the fine-tuning measure in sec. 4. Let us mention that the Markov chains are relatively short and 
therefore there is no statistical interpretation of the scatter plots 3 presented in the next sections. 

In addition to the hard cuts on the input parameter input given in Tab. 1, we restrict the 

2 Note that, in order to speed up the scan, we use the fine-tuning calculated from the tree-level potential in the 
likelihood function. After the scan, we recalculate the fine-tuning using the effective potential including one-loop 
corrections. 

3 The plots have been produced using matplotlib [19]. 
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running of the squark soft parameters by imposing 

mess 111CSS )l (3-1) 

for £ = rnQ 3 ,m^ l3 ,m^ 3 ,At,Ai ) . This avoids situations where the squark soft parameters at the 
low scale m so f t are much smaller than, for example, the loop corrections from gluinos to these 
parameters which would require an additional tuning. We do not restrict the Higgs soft parameters 
because they are already covered by the fine-tuning measure. 

We use the program NMHDECAY [21-23] contained in the package NMSSMTools 3.2.1 [24] to 
generate the sparticle spectrum. The program is modified, however, to incorporate a low messenger 
scale A mess . 

In particular, NMHDECAY calculates the soft masses Tt^, m? H and m 2 s from the input param- 
eters using the minimization conditions. Loop corrections from the Higgs-singlet sector (including 
superpartners) to the minimization conditions are neglected in this calculation. For large values 
A, k > 0(1), however, these corrections can become sizable. In order to ensure that the results 
obtained from NMHDECAY are reliable, we restrict ourselves to regions in parameter space where 
the resulting corrections to ni^j , rrijj and m| are relatively small. We therefore require that 
(cf. Eq. (2.9)) 

< |m|| . (3.2) 

The soft masses enter into our calculation of the fine-tuning measure, leading to an uncertainty of 
order 1 from the above condition. We believe that this is acceptable, given the uncertainties in the 
proper definition of the fine-tuning measure itself. In addition, NMHDECAY uses its results for 
m H u i 171 H d an d m | i n the check whether the minimum specified by the input parameters v, tan/3 
and jj, is the global minimum. We do not expect, however, that this check is significantly affected by 
the 0(l)-uncertainties in these soft masses. Nevertheless, a future implementation in NMHDECAY 
of loop corrections from the Higgs-singlet sector to the minimization conditions may be worthwhile. 
Similarly, NMHDECAY does not check whether there is a deeper minimum with VEVs v ^ 174 
GeV and vs = 0. This can be important as has recently been pointed out in Ref. [17]. 

We apply all constraints which are implemented in NMSSMTools 3.2.1, in particular several 
LEP and Tevatron exclusion limits. The program also calculates various flavor observables on 
which we impose the latest constraints. Note though that we do not try to explain the anomalous 
magnetic moment of the muon. The relic abundance of the LSP and its spin-independent proton 
and neutron cross sections are calculated using the program MicrOMEGAs 2.4.5 [23; 25]. We 
require that the universe is not overdosed and that the latest XENON100 bounds [ ] on the cross 
sections for a given relic abundance are satisfied. We take into account the constraints on oblique 
corrections and require that the Peskin-Takeuchi parameters, S and T [ '] are within their la 
limits. In addition, NMHDECA Y provides the Higgs couplings normalized to the SM as well as 
the normalized cross sections of the most important LHC Higgs detection channels. We impose 

4 In particular, via this condition the gluino limits the size of the stop mass parameters: The gluino "sucks" the 
stop mass up as has been pointed out in [20]. We find values of m.Q 3 and m U3 up to a factor 4 smaller than the gluino 
mass. This is somewhat less stringent than the ratio mg/mj < 2 given in [20], because we impose the condition (3.1) 
and include loop corrections to the stop mass parameters from m 2 Hu and A t . After imposing the LHC constraints on 
the gluino mass (see sec. 5.1), we find values of m,Q 3 and m U3 down to about 300 GeV (for sufficiently high neutralino 
mass). Such small stop mass parameters still allow for small fine-tuning. 
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dvl 



< \m 



< m 



and 



dvl 



the experimental bounds on the different Higgs cross sections. Furthermore, we restrict the scan to 
small fine-tuning by imposing T, v < 200 on the fine-tuning measure calculated from the tree-level 2 
minimization conditions. For the final results, we recalculate the fine-tuning measure including the 
dominant stop-top corrections to the effective potential. Note that we neglect the corrections from 
the Higgs-singlet sector, on the other hand, as they are restricted to be small according to (3.2). 
In order to generate data points for A mess = 100 TeV and 1000 TeV, we use the generated points 
for A mess = 20 TeV, and recalculate the constraints and fine-tuning with the new A mess . Note that 
some points are excluded for larger A mess because they have a Landau pole between 20 TeV and 
the bigger A mess . 



4 Naturalness Constraints 

4.1 Electroweak scale tuning: Large A helps 

We shall now determine to what extent the (scale-invariant) NMSSM maintains naturalness in 
light of the LHC results which place stringent lower limits on the sparticle masses. As can be seen 
from Eqs. (2.10) and (2.11), soft masses of order the TeV scale require a delicate cancellation in 
order to obtain the correct electroweak scale v = 174 GeV and thereby necessitate some amount of 
fine-tuning. 

In order to quantify this fine-tuning, we use the measure proposed in Ref. [28], which determines 
how small deviations in the input parameters of the theory affect the electroweak scale: 

dlog v 2 



= max 



L 



(4.1) 



Here the input parameters & = (m|^, , m 2 s , A, k, a\, a K , m 2 ^ , m 2 3 , m 2 ^ , A t , A^, Mi, M2, M3) are 
defined at the messenger scale A mess - This set consists of the input parameters on which the 
Higgs VEV v has a significant dependence. For the plots, we shall typically require that < 20, 
corresponding to a fine-tuning better than 5%. 

Using the minimization conditions (2.6-2.8), we can solve for the VEVs v u ,Vd and vs- We 
take the one- loop corrections to the effective potential due to third-generation quarks/squarks into 
account and choose m so ft = y/mQ 3 m U3 as the renormalization scale (cf. the discussion in sec. 2). 
The VEVs v u , and vs are then functions of running masses and couplings evaluated at this scale. 
We shall collectively denote these parameters by £i(m so ft). In order to evaluate the fine-tuning 
measure (4.1), we differentiate the minimization conditions (2.6-2.8) with respect to the £i(m so f t ) 
and solve the linearly-coupled system of equations for the derivatives du 2 /d^(m so ft). We solve the 
RG equations in the leading- log approximation to obtain an analytical expression for the £i(m so f t ) 
in terms of the parameters £«(A mess ) defined at the messenger scale. To this end, we take all relevant 
contributions at one-loop order into account and also include the two-loop corrections from gluinos 
that feed into m 2 H u , m H d an< ^ a A- Using the chain rule, we can then calculate the fine-tuning 
measure according to: 



max 



&(A mess 



V 2 d& (msoft) d& (Amess) 



(4.2) 



Let us briefly mention a subtlety when solving the RG equations: Since A and k are typically 
large and accordingly run quickly, we do not use a leading-log approximation for these couplings. 
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Instead, we solve their RG equations numerically and integrate them properly in the RG equations 
of the other couplings. Furthermore, as A and k are renormalized multiplicatively, we find that 



dln\(m soit ) 
(flnA(A mess ) 



1 , 



dlnK(m so ft) 
dln/t(A mess ) 



1 , 



dlnA(m so f t ) 
dln«;(A mess ) 







dln/t(m so f t ) 
din A(A mess ) 



. 



(4.3) 



Using these relations and the fact that the dependence of the electroweak scale on dimensionful 
parameters is already taken into account by other terms in the fine-tuning measure (4.1), we will 
approximate the logarithmic derivative with respect to A(A mess ) by that with respect to A(m so ft) 
and similarly for k. 

It is instructive to consider regions of parameter space where large stop mass parameters tUq^, 
and At at the messenger scale dominate the fine-tuning measure. In particular, the stops affect 



the Higgs mass parameter tn 2 H ^{m so ft) via the RG evolution. Using the leading-log approximation, 
we find 



m 2 Hu (m soit ) = m % u (A r _ 



[mg 3 (A mess ) + m^ 3 (A mess ) + A 2 (A mess )] log 



A r 



+ ... (4.4) 



If the contribution from the stop sector is significant, the value of m 2 H (A mess ) has thus to be tuned 
in order to obtain a value of (m so ft) not much larger than the electroweak scale. Note that 
in our scan of parameter space for A mcss = 20 TeV, we find that rrijj dominates the VEV tuning 
for approximately half of the points, while the next relevant parameters that dominate the tuning 
are A and a\. When the messenger scale is increased to A mess = 100 TeV, m# still dominates the 
VEV tuning for approximately half of the points, but the next dominating contributions are now 
M 2 and M 3 . 

We can obtain an upper bound on the stop mass parameters from naturalness by looking at how 
they affect the fine-tuning measure. Neglecting correction from the Coleman- Weinberg potential, 
we in particular find that 



d log v 2 



dlog&(A r 



3y? &(A n 



8vr 2 



log 



A r 



m so ft 



dv 2 



dm 2 



< E l 



(4.5) 



for £t = { m Q 3 T rn u 3 i A 2 ). This means that, for a given amount of fine-tuning, the parameters £t and 
thus the stop masses have an upper limit which depends on the size of the derivative dv 2 /dm 2 Hu . 
Note, however, that the logarithmic derivative with respect to m 2 H (A mess ) often gives a more 
stringent bound on stop masses from naturalness since it effectively contains the sum of the three 
individual measures in Eq. (4.5). 

The stops affect the electroweak scale v also via finite corrections which are taken into account 
by the Coleman- Weinberg potential. These have been neglected in the derivation of Eq. (4.5). 
From Eq. (4.4), it follows that as the ratio A mess /m so ft increases, the contributions from the RG 
running become more important compared to these finite corrections. For our particular choice of 
parameters, we estimate that for A mess > 100 TeV, the latter can be neglected in the calculation 
of the fine-tuning measure. 



5 This might originate from our way of generating the data points for A n 
20 TeV. 



= 100(1000) TeV compared to A n 
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Let us now discuss the derivative dv 2 /dm 2 H ^ in more detail. We first consider the situation in 
the MSSM. Neglecting loop corrections to the effective potential, we find in the limit tan /3 > 1 
that 

dv 2 _ v 2 „ / 1 



dm 2 Hu (m soit ) "m| + ° V tan /V ' ( ~ 4 ' 6 ' ) 

Thus there is no freedom in the MSSM to suppress this derivative and thereby reduce the tuning 
once the stop soft terms become large. In the NMSSM, on the other hand, we find in the limit 
A>1 that (again using the tree-level potential) 

^- r ^ V + {^). (4.7) 

The derivative is thus suppressed for large A. This effect can be observed in Fig. la, where we 
present a scatter plot of the derivative as a function of A (calculated using the tree-level potential). 
The black points have a messenger scale A mcss = 20TeV, while the orange (yellow) points have 
Amess = 100 (1000) TeV. There are fewer points for larger A mess because some points are excluded 
by a Landau pole below A mess . The suppression at large A allows for larger stop masses compared 
to the MSSM for a given amount of fine-tuning. In Fig. lb, we show a scatter plot of the stop mass 
mj versus A, where all points have a fine-tuning better than 5%. Indeed, the maximal value of 
that we find grows with A, corresponding to the suppression of the derivative dv 2 /dm 2 H with A. In 
particular, naturalness can be maintained for stop masses m t - > 1 TeV as opposed to the MSSM 
where lighter stops mi < 600 GeV are required. Note that when k ~ A, the suppression in (4.7) 
becomes 1/A 2 , and agrees with the result in [29]. 

In addition, we find that the stop/top-contributions to the Coleman- Weinberg potential alleviate 
the fine-tuning. This has been noticed before in the context of the MSSM [ .]. Indeed, these 
corrections increase the quartic Higgs coupling, resulting in a reduced sensitivity of the electroweak 
scale in analogy to the case of large A. We find that on average, this effect reduces the fine-tuning 
measure by 10% — 20% compared to the case when only the tree-level potential is taken into account. 

This observation, however, implies a caveat: We have neglected the contribution from the Higgs- 
singlet sector to the Coleman- Weinberg potential in the calculation of the VEV fine-tuning measure. 
In the region of large A, this contribution is potentially important and could modify the suppression 
of the derivative dv 2 / dm 2 H ^ with A. This is indeed what we expect in the decoupling (or SM) limit 
when the Higgs potential can be described by a one-Higgs-doublet model: V ~ — m\ eS h 2 + \ e gh 4 . 
Here m\ eff is a combination of the dimensionful parameters in the potential, whereas A e g is deter- 
mined by the gauge couplings and A. Both m? eS and A e ff also receive important loop corrections. 
Minimizing this potential, we find that the Higgs VEV is given by 

2 

2 m h,cS Q , 

v ~ — - J — . (4.8) 

A e ff 

From this relation, we again see that large A can alleviate the fine-tuning problem. Indeed, in- 
creasing A typically also gives a larger effective quartic coupling A c fr • This in turn means that the 
dimensionful parameter m\ cS is larger and can accordingly also tolerate larger loop corrections 
without affecting the Higgs VEV too much. In the decoupling limit, however, the dimensionful 
parameter in the potential is essentially the Higgs mass (up to an O(l)-factor). With a fixed Higgs 
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(a) (b) 

Figure 1: Scatter plots of (a) the tree-level derivative \dv 2 /dm 2 H J\, and (b) the lightest stop mass 
wif,, as a function of the Higgs-singlet coupling A. In (a) the black, orange, yellow points correspond 
to A 

m ess — 20, 100, 1000 TeV, respectively. All points in (b) have A mess — 20TeV and a tuning in 
the Higgs VEV better than 5%. In (b) the green points have a combined tuning (cf. sec. better 
than 5%, i.e. Yj h Y? < 20, for the blue points it is between 1% and 5%, while for the red points it is 
worse than 1%. The derivative dv 2 /dm 2 H is suppressed for larger values of X, allowing for mi to 
become as large as 2.5 TeV for a combined tuning better than 1%. All points satisfy the constraints 
discussed in sec. 5. 

mass (and thus m^eg ~ 126 GeV), a larger coupling A does therefore not allow a larger dimen- 
sionful parameter m/j, e ff in the potential and should accordingly not alleviate the fine-tuning in the 
decoupling limit. 

We believe that loop corrections from the Higgs-singlet sector may play an important role 
in this context. Indeed, once the VEV and mass in Eq. (4.8) are fixed, the effective quartic 
coupling A e ff is fixed as well. This means that, at large A, an accidental cancellation between 
the A-contribution to A c ff and the loop corrections has to occur in order to bring A e g down to 
the required value. We will discuss this tuning (which can be phrased as a tuning in the Higgs 
mass) in some detail in the next section. Given that corrections from the (s)top sector raise the 
quartic coupling (or, equivalently, raise the Higgs mass), these corrections can not be responsible 
for this cancellation. The contribution from the Higgs-singlet sector, on the other hand, can lower 
the quartic coupling (or, equivalently, lower the Higgs mass) [2]. We therefore expect that these 
corrections may counteract the suppression of the derivative dv 2 / dm 2 Hu at large A in the decoupling 
limit. Let us emphasize, however, that most of our points clearly deviate from this limit, where the 
potential and thus the effect of large A is more complicated. In addition, the presence of additional 
(s)particles with masses 0(v) can lead to important non-trivial VEV-dependent contributions from 
the Coleman- Weinberg potential to Eq. (4.8) even when the Higgs couples very SM-like. In any 
case, it may be worthwhile in the future to include the contributions from the Higgs-singlet sector 
to the Coleman- Weinberg potential in the calculation of the fine-tuning measure. 
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4.2 Higgs mass tuning: Large A hurts 

In the preceding section, we have discussed how a large Higgs-singlet coupling A relieves the tuning 
which is required in order to obtain a Higgs VEV at the electroweak scale. As we will now see, 
however, large A introduces a new source of tuning associated with obtaining a Higgs mass around 
126 GeV. 

Let us first consider the limit in which the Higgs (i.e. the state that we identify with the newly 
discovered resonance at 126 GeV) couples exactly like the SM Higgs. As discussed in sec. 2, this 
is the case if the Higgs corresponds to the linear combination h! in the basis (2.18). Its mass at 
tree-level is then given by (2.19). There are typically important loop corrections to this tree-level 
mass which come dominantly from the top-stop and Higgs-singlet (including superpartners) sector. 
We shall denote these contributions respectively by 5rnf lstop and Smf lS . Neglecting subdominant 
loop corrections, the Higgs mass can be written as 

m 2 h ~ m| cos 2 2/3 + AV sin 2 2/3 + 6m 2 h>stop + 5m 2 h S . (4.9) 

As we will discuss in sec. 5.2, electroweak precision tests require small tan/3 if A is relatively large. 
The tree-level mass thus goes like m/j )tree ~ Xv in this regime. In particular, it becomes too large 
for A > 1. Large loop corrections are then required to bring the Higgs mass down to 126 GeV. As 
is well known, the stop-top sector raises the Higgs mass and is thus not suitable for this purpose. 
The Higgs-singlet sector, on the other hand, can give a sizable negative contribution [2]. 

In our scan, we find values of A up to 2.4 and tan/3 close to 1. The tree-level mass-squared 
for such a large A and small tan/3 becomes Tnf li ~ 10 x (126 GeV) 2 , i.e. a factor 10 larger than 
what is required. In order to obtain the correct Higgs mass, a cancellation at the 10%-level has 
to occur between this part and the loop corrections. The input parameters of the theory have to 
be tuned accordingly. Note that this type of tuning is different from (and thus comes in addition 
to) the tuning which is required to obtain the correct Higgs VEV. Indeed, the tree-level mass in 
Eq. (4.9) is expressed in terms of the Higgs VEV (and tan/3). The tuning in the Higgs VEV only 
ensures that v = 174 GeV in Eq. (4.9). The hierarchy for large A between m^tree with this v and 
the measured value of 126 GeV necessitates an additional tuning in the Higgs mass. 

The discussion so far has focused on a Higgs that couples exactly as in the SM. The range of 
Higgs signal strengths which are consistent with the measurements at ATLAS and CMS, however, 
still allow for some admixture of the orthogonal linear combination H' and the singlet s to the 
Higgs. In our scan, we find that the state which we identify with the resonance at 126 GeV always 
coincides with the lightest CP-even Higgs. The admixture of H' and s then gives an additional 
negative contribution to eq. (4.9) due to the usual effect of level repulsion. If the admixture of 
the orthogonal doublet is negligible, the contribution to the Higgs mass due to the singlet can be 
approximated as [17; 31] 



sin 2/3 



1 + 



2Xkvs 



Sm 2 ^ ~ -AV ax —^: 

1 + 1^ 1 8K 2 tr| 



(4.10) 



This expression becomes exact in the limit ml/h'/M^ ~~ ^ 0j wnere ^33 denotes the (3, 3)-element 
of the CP-even mass matrix given in (2.14) and myh' is the tree-level mass in the SM-limit given 
in (2.19). We find that the approximation m\ = m 2 , ft , + 5m 2 mix with the above mixing term 
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reproduces the Higgs mass at tree-level (obtained from diagonalizing the CP-even mass matrix in 
(2.14)) to better than 10% for the vast majority of points in our scan. This corresponds to the 
fact that the admixture with the orthogonal doublet and thus its contribution to the Higgs mass 
is typically small in our sample of points. 

The 'pull-down' of the Higgs mass due to the admixture with the orthogonal doublet and the 
singlet was considered in [17] as a means to obtain the correct value of 126 GeV in the regime 
of large A. We find, however, that this admixture modifies the mass by at most 40% relative to 
(2.19) for almost all our points. The 'pull-down' at tree- level is thus typically not a large effect. 
Instead, loop corrections tend to dominate and we find that Sm^ stop + <5m| g > 2 x Sm^ mix for 
almost all our points. We believe that this is due to the following reasons: Firstly, in the regime of 
large A, loop corrections from the Higgs-singlet sector are important (the latter were neglected in 
[17]). Secondly, the Higgs signal strengths which were measured at ATLAS and CMS already put 
a significant limit on deviations from a SM Higgs. This restricts the admixture of H' and s to the 



Higgs and thus the size of 5m\ 

We thus find that, at least for large A, some accidental cancellation among different contributions 
to the Higgs mass is necessary in order to obtain the correct mass. It is useful to quantify this 
tuning in the Higgs mass in analogy to the fine-tuning measure for the Higgs VEV. We therefore 
define 



mix' 



max 



d log m| 



dlog& 



(4.11) 



where £j = (A, n, ax, a K , TriQ 3 , m 2 3 , rn^ 3 , At, A^, Mi, Mi, M3) are input parameters that affect the 
Higgs mass. Note that we do not include the soft masses m H d anc ^ m s among these parame- 

ters. Instead, we consider a basis in which these soft masses have been replaced by the Higgs VEVs 
v u , Vd and vs- We then keep the latter fixed when evaluating (4.11). This allows us to separate 
the fine-tuning in the Higgs mass from that in the Higgs VEV. As before, a measure T, h = 20 
corresponds to a tuning of 5%. 

We present a scatter plot of Y* h as a function of A in Fig. 2a. 6 Notice that the smallest tuning 
in the Higgs mass that we find for a given A grows approximately like A 2 . This can be understood 
from the fact that the tree-level part of the approximate mass formula (4.9) grows like A 2 . A 
correspondingly larger cancellation among this tree-level part and the loop correction is therefore 
required at larger A in order to bring the Higgs mass down to 126 GeV. We find that less tuning 
is needed, on the other hand, for larger tan ft. This can be seen in Fig. 2b and is due to two 
effects: The tree-level mass becomes smaller for larger tan f3 if A is kept fixed (assuming that A 
is so large that the term m 2 z cos 2 2/3 in (4.9) is subdominant). In addition, larger values of tan/3 
require smaller values of A in order to satisfy constraints on the T-parameter (see the discussion in 



In order to evaluate Eq. (4.11), we calculate the tree-level contribution m h trcc to the Higgs mass from the CP-even 
mass matrix in (2.14) and the loop correction <5m| loop using NMSSMTools. We then make the approximation 

dlogm 2 h & drnl itlee 



d log £i m 2 h d& 

i.e. we neglect the derivative of the loop correction 8m\ loop with respect to the parameters £j and thus a possible 
tuning among those & which enter into 5?n^ loop . This approach instead only takes into account a tuning among 
parameters within m| tree and an accidental cancellation between 171% tree and 8m\ loop . Given that the tuning in 
the Higgs mass is typically of the latter type as discussed around Eq. (4.9), we believe this to be a reasonable 
approximation. 
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Figure 2: The Higgs mass tuning measure T, h as a function of (a) A and (b) tan/3 where A mess = 
20TeV and the Higgs VEV tuning is better than 5%. The green, blue and red points correspond to 
a combined tuning better than 5%, between 1% and 5%, and worse than 1%, respectively. Notice 
that the amount of tuning grows approximately like A 2 in (a) and that it is reduced for larger tan /3 
in (b). All points satisfy the constraints discussed in sec. 5. 



sec. 5.2). This gives another suppression on the tree-level part of the Higgs mass. 

We thus find that the Higgs mass often is a significant additional source of fine-tuning in 
the scale-invariant NMSSM. This type of tuning should accordingly be taken into account when 
assessing the naturalness of the model. In the following, we will adopt a simple possibility and 
multiply the tuning measures for the Higgs VEV and the Higgs mass to quantify the combined 
tuning. ' This approach is justified since T, h is constructed in a way that makes it independent of 
the tuning in the Higgs VEV as discussed before. 

In Fig. 3a, we show a scatter plot of the combined tuning as a function of A. In addition to the 
messenger scale A mess = 20TeV mainly considered in this paper (black points), we have included 
points with A mess = 100 TeV (orange) and A mess — 1000 TeV (yellow). Note that the combined 
tuning increases with growing messenger scale. This is expected since loop corrections to soft 
parameters from RG running also increase with the messenger scale (cf. Eq. (4.1)). Furthermore, 
notice that there is a minimum for the tuning at values A ~ 1. This can be understood as follows: 
Let us first consider the region A > 1. We expect that the smallest tuning measure which is 
achievable for a given A, decreases with increasing A in this region. This is indeed what we observe 
in our sample of points. The trend, however, is less pronounced than the A _2 -dependence that we 
may naively expect from the arguments in the preceding section. 8 Since T, h grows approximately 
like A 2 , on the other hand, this explains the growth with A of the combined tuning in the region 
A > 1 of the scatter plot. In the region A < 1, the tuning measure T, h becomes approximately 
independent of A as can be seen in Fig. 2a. This corresponds to the fact that the Higgs mass 
at tree-level is below 126 GeV for A < 1 so that no tuning arises in the Higgs mass. Since loop 

7 In a statistical sense, if two quantities are not correlated, the probability involving both is P(A(lB) — P(A)*P(B). 
8 We note that this partly results from our linear scanning of the parameter space (see sec. 3) which prefers larger 
values for the soft masses and thereby counteracts the suppression of E" for large A. 
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Figure 3: The combined tuning S^S'' as a function of (a) the Higgs-singlet coupling A and (b) m so f t . 
The points in (b) have A mess = 20TeV and a tuning in the Higgs VEV better than 5%. The green, 
blue and red points in (b) correspond to a combined tuning better than 5%, between 1% and 5%, 
and worse than 1%, respectively. The brown band in (b) corresponds to the Higgs VEV tuning for a 
generic region of the MSSM with G [124, 127] GeV. All points satisfy the constraints discussed 
in sec. 5. 



corrections from the top-stop sector are then required to raise the mass, the tuning measure H v 
in turn grows for decreasing A. Together this explains the behavior of the combined tuning in the 
region A < 1 of the scatter plot. 

In Fig. 3b, we show a scatter plot of the combined tuning as a function of m so f t = ^/m.Q 3 m U3 . The 
brown band maps the tuning in the MSSM, where we have chosen A mess = 20 TeV, /i = 200 GeV, 
tan /3 = 20 and m a = 1 TeV for the mass of the CP-odd scalar to represent a generic region of 
parameter space. Furthermore, the A-term At is varied such that nih = 124 (127) GeV for the lower 
(upper) boundary. We have used the program FeynHiggs [32] to calculate the Higgs mass. Note 
that there is no tuning in the Higgs mass required in the MSSM. Indeed, the mass at tree-level 
is always too small and large loop corrections are required to lift the mass to 126 GeV. We have 
accordingly set S ft = 1 for the MSSM. We see that, even taking the combined tuning into account, 
there is a large region of parameter space for which the (scale- invariant) NMSSM requires less 
tuning than the MSSM. This is not surprising as values A ~ 1 do not introduce a significant tuning 
in the Higgs mass but alleviate the need for heavy stops to raise the mass as in the MSSM. 

There can also be a tuning associated with the lightest CP-odd Higgs state. This can be seen 
by noting that the lightest mass-squared eigenvalue in (A.l) is given by 

2 Xkv 2 . . 

m a, - -3a K v s + 9 — iy-axvs ■ (4.12) 
m A 

Clearly if m ai < 100 GeV with soft mass parameters in the TeV range, a tuning will occur, which 
was verified numerically: a tuning in m ai at the 5% level is found for m ai < 100 GeV for parts 
of the parameter space. We also checked numerically that all of the other particle and sparticle 
masses in the regions of interest have tunings better than 50%. So we conclude that if at some 
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time in the future a CP-odd scalar is discovered with mass m ai < 100 GeV it could result in an 
additional source of tuning from the perspective of the model under consideration. 



5 Limits from Colliders and Cosmology 

Physics beyond the SM is already strongly constrained by direct searches at colliders, as well as 
indirectly by electroweak precision tests and flavor physics. Several such limits (in particular from 
LEP, the Tevatron and the B-factories) are already implemented in NMSSMTools 3.2.1. We apply 
all these limits in our scan, except for the flavor constraints for which we use updated results. 
Also note that we do not attempt to explain the discrepancy in the measurement of the anomalous 
magnetic moment of the muon, [g — 2)„, with the SM expectation. Furthermore, we require that 
the LSP abundance does not overdose the universe and that it satisfies direct detection limits. 

In addition to the limits implemented in NMSSMTools 3.2.1, we take into account several 
searches at the LHC as well as electroweak precision tests. In the following, we focus on our 
additional constraints and refer the interested reader to the documentation of NMSSMTools 3.2.1 
for the remaining limits. 



5.1 Direct limits from colliders 

We identify the bosonic resonance, which was recently discovered in the diphoton and diboson 
channels at ATLAS [ !] and CMS [3 ], with one of the three CP-even states in the Higgs sector of 
the NMSSM. We denote this state by h. In agreement with the current experimental uncertainties, 
we require that its mass lies in the range 124 GeV < < 127 GeV. Furthermore, we require that 
its signal strength is within the la range measured at the LHC for the most constraining channels. 
More precisely, we take the average value of the signal strengths quoted by ATLAS [33; 35; 36] and 
CMS [34; 37] for each channel and combine the errors in quadrature. Here the signal strength is 
defined as 

_ a(h) x BR(h -> X) 
X a(h SM ) x BR(h SM -> X) ' y ■ > 

where hsu denotes the SM Higgs, a the production cross section at the LHC and BR the branching 
fraction to final states X. We then require that 

0.81 < R ZZ < 1-32, 0.74 < R ww < 1.40, 

0<R bi < 1.10, 0.27 < R TT < 1.15, (5.2) 

where for R bb , the Higgs production is assumed to be in association with a vector boson. However, 
in order to allow for a Higgs with SM properties in our scan, we do not impose the corresponding 
range for the diphoton channel: 

1.40 < Ry 7 < 1.99 . (5.3) 

Instead, we discuss it separately in sec. 6. 

The Higgs searches by ATLAS [33] and CMS [34] also constrain the heavier CP-even states S2 
and S3. At high masses, decays into WW and ZZ dominate. These channels are highly constrained 
and the production cross section times branching fraction (into WW and ZZ) of these states has 
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to be significantly smaller than the corresponding rate of a SM Higgs at the same mass [38; 39]. 
More precisely, we shall require that 



a(si) x BR( Si -> ZZ) 



< 0.09 



a( Sl ) x BR{ Si -»■ WW) 



< 0.2. 



(5.4) 



a(h SM ) x BR(h su -> ZZ) 



ct(/ism) x 5fl(/i SM -> WW) 



Besides limits on CP-even Higgs bosons, we impose all bounds on sparticle masses coming from 
LEP and the Tevatron. In particular, these bounds are (cf. [40]) > 95.7 GeV for the stop, 
m~ b > 89 GeV for the sbottom, m H ± > 79.3 GeV for a charged Higgs, and m^± > 94 GeV for a 
chargino. 

Furthermore, we take into account the latest SUSY searches at the LHC. The ATLAS and CMS 
experiments have presented exclusion plots for sparticle masses in various simplified models, where 
most of the sparticles are decoupled. Note that, for a typical point in our parameter space, more 
sparticles will be light than assumed in these simplified models. This can lead to longer decay 
chains as those considered in these studies and could weaken certain limits or make them more 
stringent. In our scan, we shall nevertheless impose bounds derived from these simplified models. 
We do not expect that our results are affected by this simplified treatment of the LHC limits for 
the following reasons: Firstly, we impose bounds which are more conservative than those coming 
from the simplified models (the precise limits will be given below). This leaves some room in case 
the simplified-model limits are not sufficient to ensure that a given point in our parameter space 
does not lead to a visible excess in the LHC searches. Secondly, we are mostly interested in the 
upper limits on sparticle masses coming from naturalness, whereas direct searches give lower limits. 
Therefore, we expect that our results do not depend sensitively on the latter. 

The experimental limits depend on the nature of the LSP. Let us therefore first discuss the 
different possibilities for the LSP in our scenario. In order to minimize the size of loop corrections 
to soft masses, we require that soft breaking terms are generated at a low scale A mess (which we refer 
to as the 'messenger scale') . One can imagine various types of UV completions above this scale. 
For example, A mess could be associated with the mass of messenger fields like in gauge-mediated 
SUSY breaking. In this case, a gravitino with mass m 3 / 2 ~ w>softA m ess/A^Pl ~ 0.01 eV would be the 
LSP. Alternatively, A mess could be a cutoff at which the Higgs sector and/or other NMSSM fields 
emerge as composites of an underlying strongly coupled theory. This is, for example, realized in 
the models of [7-10]. This case allows for the possibility that the LSP resides in the visible sector 
instead of being the gravitino. Early-universe cosmology then favors a neutralino as the LSP. 

In our scan, we shall focus on the case of a neutralino LSP. This is partly motivated by our 
desire to study the prospects of dark matter in this scenario (gravitinos with mass 7713/2 ~ 0.01 eV 
are underabundant). In addition, there are currently only a few searches which consider a spectrum 
with a gravitino LSP (usually in the framework of general gauge mediation) and decoupled first- 
two-generation sparticles. This makes the implementation of LHC constraints without performing 
detailed simulations difficult. 

We believe, however, that our results are relatively independent of this choice of the LSP. This 
is for the following reasons: We have performed a dedicated scan in which we made no requirements 
on the nature of the lightest superpartner (i.e. the NLSP if the gravitino is the LSP). After imposing 
the limits discussed so far in this section and those from LEP and the Tevatron, we have found 
that in 98% of the cases the LSP is a neutralino (with various admixtures of wino, bino, singlino 
and Higgsino components). Requiring that the LSP is a neutralino, as we do in our scan, thus 
does not eliminate many points. There are dedicated searches for gluinos and stops in scenarios 
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with a gravitino LSP, decoupled first-two-generation sparticles and a bino-like neutralino NLSP 
[41-43] (in events with diphotons plus $t) or a Higgsino-like neutralino NLSP [42] (in events with 
a Z boson, jets plus fir)- The limits on stops and gluinos masses are comparable to those for a 
neutralino LSP except that the constraints are almost independent of the mass of the neutralino 
NLSP (in particular nig < 1.3 TeV is excluded for all NLSP masses m^o < nig). We therefore 
expect that the set of points which satisfy the current LHC limits does not differ significantly for 
the cases of a gravitino and neutralino LSP. Finally, dark matter constraints eliminate some points 
for a neutralino LSP (due to overclosure and direct detection limits as discussed in sec. 5.4). No 
such constraints apply for a gravitino LSP with mass m 3 /2 ~ 0.01 eV as its abundance is sufficiently 
small. However, we have found that less than 1% of the points which survive all other limits are 
eliminated by the dark matter constraints in the former case. For these reasons, we expect that the 
set of points which pass all constraints would not change significantly if we considered a gravitino 
LSP instead of a neutralino LSP. 

Let us now discuss the LHC limits that we impose for the case of a neutralino LSP. An important 
search by ATLAS uses final states with b-jets and missing transverse energy {fir) in 12.8 fb _1 of 
data taken at 8 TeV [44]. Bounds are derived in two simplified models where gluinos are directly 
produced and decay to the neutralino LSP and either bb or ti (with 100% branching fraction). 
All sparticles except the gluino and the neutralino are decoupled. Since the limits from the decay 
g bbxi are slightly more stringent than those from g — > tixi, we will impose the former. Following 
the exclusion curve presented in [44] for this case, we shall conservatively exclude the region 

m g < 1310 GeV if m^o < 650 GeV. (5.5) 

Another search by ATLAS in events with jets, leptons and $t [45] gives weaker constraints on the 
gluino mass. Furthermore, CMS has used less data for its gluino searches than ATLAS and the 
resulting limits are correspondingly less stringent too. 

An important search by CMS uses final states with b-jets and J^t in 11-7 fb^ 1 of data at 
8 TeV [46]. Constraints are derived on a simplified model, where the sbottom is directly produced 
and decays with 100% branching fraction into the neutralino LSP and a b-quark. Again all other 
sparticles are decoupled. Following the relevant exclusion curve presented in [16; 47], we will forbid 
points with 

150 GeV < m~ b < 650 GeV if m^o < 230 GeV . (5.6) 

Note that ATLAS did not update their sbottom search using 8 TeV data yet and their bounds 
are correspondingly weaker. An analogous simplified model with stops instead of sbottoms is 
constrained by two ATLAS searches: In [48], final states with jets and JGt are analyzed using 4.7 
fb _1 of data at 7 TeV. In [ ], these final states are required to also have an isolated lepton. 
Following the exclusion plots in these papers, we shall conservatively forbid points in the region 

220 GeV < m~ t < 500 GeV if m^o < 160 GeV . (5.7) 

Note that, even though the latest stop search by CMS [50] uses 8 TeV data and higher integrated 
luminosity, it results in comparable limits. 

Final states with leptons and $t are used in dedicated searches for neutralinos and charginos 
by ATLAS [51] and CMS [52; 53]. Limits are derived in a simplified model where the chargino- 
neutralino pair X1X2 1S produced and decays with 100% branching fraction to Xi^xlW^ . Further- 
more, to simplify the parameter space, it is assumed that the LSP x.i is bino-like, xfand x.2 are 
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wino-like and that the latter two particles have the same mass. The more stringent limit comes 
from the CMS search [53] which excludes approximately the region 

m-± = m^o < 330 GeV if ra^o < 120 GeV and m^o + mw < m~± . (5.8) 

Xl X2 Al Al Xl s ' 

The limit does not apply for m^o + mw > m-± because on-shell Ws and Zs are assumed in the 

Xi Xi 

CMS search. We find that this does not pose any constraint on our generated data points, because 
there are only a few points with a dominantly bino LSP and xfiX® being wino-like. All of these 
points have an LSP mass m^o > 120 GeV. 

5.2 Electroweak Precision Tests 

The additional particles of the NMSSM compared to the SM contribute to the Peskin-Takeuchi 
parameters [27] S and T. Using the latest measurements, the Particle Data Group [40] quotes 
the ranges S = -0.04 ± 0.09 and T = 0.07 ± 0.08 at 95% C.L. for a reference Higgs mass of 
nT>h,ref = 117 GeV with a correlation between S and T of 88%. We constrain the model parameters 
by demanding accordance with these measurements at 95% C.L. . 

We can distinguish three sectors that contribute to S and T: i) Higgs and singlet scalars, ii) 
stops and sbottoms and iii) neutralinos and charginos. Since the expressions for the contributions 
of these sectors can be found in the literature, we shall not repeat them here. We use the formulas 
given in [54-56] for the neutralino-chargino and stop-sbottom sector and those from [3; 12] for the 
singlet and Higgs scalars. In order to compare with the values quoted by the Particle Data Group, 
the contribution of a SM Higgs with reference mass mh, r ef = 117 GeV has to be subtracted from 
the latter. 

Note that we do not decouple the electroweak gauginos but take all charginos and neutralinos 
into account, in contrast to the analysis in [3; 12]. In these references, it was found that accordance 
with electroweak precision tests requires tan/3 < 3 for large A ~ 2. Here the largest contributions 
come from the stop-sbottom and chargino-neutralino sector. We have verified that we recover a 
similar constraint on tan/3 in the limit of decoupled electroweak gauginos and A ~ 2. Smaller 
A allows for larger tan/3 and in our scan, we find values up to tan/3 ~ 4.2 for A ~ 1.1. This 
dependence on A is related to the breaking of the custodial SU(2) for tan /3 ^ 1 (which in turn 
leads to a correction to the T parameter). Indeed, the Higgsino-singlino sector can provide sizable 
contributions to the T parameter due to large off-diagonal terms in the mass matrix (see [12] for a 
discussion of this in the limit of decoupled electroweak gauginos) . This mixing and thus the effect 
of custodial symmetry breaking due to tan/3 > 1 is enhanced for large A. 

5.3 Flavor Constraints 

The package NMSSMTools allows to calculate several flavor observables. We use the calculated 
values to impose the current constraints. For the mass differences in the B° system, AMj, and the 
Bg system, AM S , we use the current 2a ranges from the HFAG group [57]: 

A M s = (17.719 ± 0.086) ps _1 , 
AM d = (0.507 ± 0.008) ps" 1 . 
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We also impose the 2a ranges on the branching ratios of the B decays B + — > t + v t and B —> X s ^f 
from the same group: 

Br(B + -> t + u t ) = (1.67 ± 0.60) x 10 -4 , 
Br(5 -> AT s7 ) = (3.55 ± 0.48 ± 0.18) x 10" 4 . 

In case of the rare decay B® — > we use the 2a ranges from the recent LHCb measurement [58]: 

Br(B° -+ = 3.7**£™ x 10~ 9 . 

Since small values of tan/3 ~ 1 are needed in order to have a sizable tree-level contribution 
to the Higgs mass and to satisfy electroweak precision tests in our model, the dominant flavor 
constraints arise from the contribution of the charged Higgs to B — > X s ^f. Indeed, this contribution 
dominates in the small tan /3 limit when \iA% > 0, while other flavor changing effects become only 
relevant when tan/3 S> 1 [59]. 

We do not attempt to explain the discrepancy in the measurement of the anomalous magnetic 
moment of the muon, (g — 2)^, with the SM expectation. Furthermore, we do not consider any 
constraints from lepton flavor violating processes since we assume decoupled sleptons, mr « A mess . 

5.4 Cosmological Constraints and Dark Matter 

As discussed in sec. 5.1, we focus on the case of a neutralino LSP in our scan. Since we assume 
i?-parity conservation, we have to ensure that the LSP relic abundance does not overdose the 
universe. Using the WMAP-7 measurements of the dark matter abundance [60], we require that 
ft LSP /i 2 < OwMAP-7/1 2 = 0.1120 ± 0.0056. In addition, we impose the latest limits on direct 
detection from XENON100 [ ]. We calculate the relic abundance and the direct-detection cross 
sections using the program MicrOMEGAs [23; 25]. 

Note, however, that constraints from overclosure and direct detection are avoided if the gravitino 
is the LSP. Indeed, if the scale A mess is associated with the mass of messenger fields like in gauge 
mediation, we expect the gravitino to have a mass m 3 / 2 ~ 0.01 eV (see sec. 5.1). The goldstino 
component of the gravitino couples with strength l/(m 3 / 2 Mp), where Mp is the reduced Planck 
mass. Such a light gravitino is therefore in thermal equilibrium with the plasma in the early 
universe [61]. After freeze-out, its relic abundance is f2 3 / 2 /i 2 « (m 3 / 2 /keV) x (100/<?*(7/)), where 
g*(Tf) is the effective number of degrees of freedom at freeze-out (typically g* ~ 100 — 200) [62]. 
For m 3 / 2 ~ 0.01 eV, we thus find that $7 3 / 2 /i 2 <C $7wmap-7^ 2 - This also avoids all constraints 
from direct detection. Furthermore, we note that the available collider limits for a scenario with 
neutralino NLSP and gravitino LSP are not significantly different from the limits for a neutralino 
LSP (see sec. 5.1). In this section, we shall therefore include points which would be excluded for 
a neutralino LSP and will interpret them in the context of a scenario with neutralino NLSP and 
gravitino LSP. 9 In particular, we require the gluino for these points to satisfy the limit m g > 1.3 TeV 
(cf. sec. 5). 

The scatter plot in Fig. 4 shows the relic density Qh 2 of the neutralino LSP as a function of 
its mass. We have colored the points according to the species which dominates the neutralino 
composition: green, blue, orange, and red points correspond to neutralinos which are respectively 
mostly singlino, Higgsino, wino or bino. Purple points, on the other hand, are excluded due to 

Alternatively, these points could be included if the reheating temperature is sufficiently low (see e.g. [63]). 
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Figure 4: A log-linear plot of the LSP relic density Qh 2 versus the LSP mass Tn^o. The colors 
denote the dominating component of the neutralino: green, blue, orange, red corresponds to a lightest 
neutralino which is dominantly singlino, Higgsino, wino or bino, respectively. The purple points are 
excluded by the WMAP-7 or XENON100 constraints, but they are still viable for a scenario with 
gravitino LSP and neutralino NLSP. The gray band depicts the experimentally measured DM relic 
density. All points satisfy the constraints discussed in sects. 5.1 to 5.3. 

overclosure or direct detection for a neutralino LSP and accordingly require a gravitino LSP (in 
which case the neutralino NLSPs decay on short time scales to gravitinos and SM particles). It 
turns out that the neutralino NLSP for these points is mostly bino. 

Note that for most of our points, the neutralino LSP is underproduced compared to the relic 
density which is required to account for the dark matter. This meets our expectations for a 
neutralino which is mostly wino or Higgsino, since its couplings to the Z-boson and the Higgs 
then yield a large annihilation cross section. It was found in [3; 12], on the other hand, that the 
right relic abundance can be obtained for a mostly-singlino neutralino. It is possible that our 
scan missed such points because of the scanning procedure or the fine-tuning constraint. We also 
note that ref. [3] considered a superpotential which is at most quadratic in the singlet superfield S 
(ref. [12] assumed the scale-invariant NMSSM though). The presence of the cubic term in the scale- 
invariant NMSSM allows for new annihilation channels and results thus in a smaller relic abundance 
for a mostly-singlino neutralino. We find that the CP-even scalar S2 (the second-lightest after the 
Higgs) is typically mostly singlet and not very heavy, m S2 < 1.1 TeV. The cubic interaction in 
the superpotential then allows for strong s-channel annihilation of a mostly-singlino neutralino via 
S2- Similarly, if the Higgs has a significant singlet admixture, the Higgs can mediate the s-channel 
annihilation. 

We find a few points for which the right relic density is obtained to account for the dark matter. 
For these points, the neutralino LSP is mostly bino with a small admixture of wino and Higgsino, 
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(a) (b) (c) 

Figure 5: Main annihilation channels of the neutralino LSP for points with relic density 
Q^o h 2 ~ 0.1. Though X\ * s mostly bino for these points, its Higgsino component provides the 
strongest annihilation channels. 

whereas the singlino component is negligible. Its mass lies in the range < m^o < 500 GeV 

(the lower limit results from the constraint on the invisible decay width of the Z). Furthermore, we 
find that the lightest chargino is mostly Higgsino with mass m~± > 200 GeV, whereas the charged 
Higgs has a mass m H ± > 750 GeV. We therefore expect that t-channel annihilation of the LSP 
via a chargino and LSP-chargino co-annihilation is suppressed due to the heavy charged Higgs. 
The second-lightest neutralino, on the other hand, is mostly singlino with a small wino component 
and a mass m^,o > 350 GeV. This makes significant co-annihilation of this state with the LSP 
unlikely. Stops and sbottoms also have large masses, g > 1 TeV, which suppresses annihilation 
channels involving squarks. Instead, we find that the most relevant annihilation channels result 
from the small but non-vanishing Higgsino component of the neutralino LSP. This allows for s- 
channel annihilation via a\ to ha\ or via the Higgs to 66, WW, where the Higgs is typically on 
resonance, 2m^o ~ nih- The corresponding Feynman diagrams are shown in Fig. 5. The fact 
that the Higgsino component of the neutralino LSP is rather small is compensated by the large 
A-coupling. To recapitulate, we find that the right relic density can be obtained for a neutralino 
LSP if it is mostly bino and has a small but non-vanishing Higgsino component. 

6 Phenomenological Implications 

6.1 Particle spectrum 

The spectrum of particles in the scale-invariant NMSSM that plays a dominant role in the Higgs 
sector is constrained by two competing effects. Naturalness prefers the Higgsino, stop/sbottom 
and gluino to be light, while experimental constraints, mostly from the 8 TeV LHC run, place 
increasingly sizable lower limits on these sparticle masses. Nevertheless there still remains a range 
of particle masses that has not yet been fully explored at the LHC where the tuning is not too 
severe. 

In Fig. 6, we show the ranges of masses in our scan that satisfy all constraints described in 
sections 3, 4 and 5. An effective /i-term in the range fx G [260, 840] GeV is generated. The state h 
in the mass window [124,127] GeV (that we identify with the newly discovered bosonic resonance 
at the LHC) always corresponds to the lightest CP-even state. If we demand a total tuning 
better than 1%, this state is mostly h u in its composition. Notice that usually independently of 
tuning, most particles in the Higgs-singlet sector and part of the neutralino-chargino sector can 
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Figure 6: The particle spectrum of the scale-invariant NMSSM with A mess = 20 TeV that we find in 
our scan satisfying all constraints described in the text. The green, blue and red points correspond 
to a combined tuning (T> h Yi" ) better than 5%, between 1% and 5%, and worse than 1%, respectively. 
All points satisfy the constraints discussed in sec. 5. Note that there are correlations among the 
masses which are not shown in the figure, but are described in sec. 5. 

have masses below 1 TeV. Most prominently, we find such low masses for the states S2, a\, ai and 
X?) x!L xf> respectively. In particular, the lightest chargino for a combined tuning better than 1% 
has a mass m~± < 575 GeV which can lead to loop-induced enhancement in the Higgs diphoton 
decay (cf. sec. 6.2). The colored sector tends to be heavier but the possibility of stops, sbottoms 
and gluinos in the windows £ [220, 1000] GeV, mg £ [216, 1000] GeV, nig > 1.2 TeV remains 
viable. In contrast to the MSSM however, the lightest stop is not necessarily light. We find stop 
masses m£ up to 2.5 TeV for a combined tuning better than 1%. Given that in the absence of new 
colored states, the colorless states are only produced via electroweak processes, it will therefore 
require searches at the 14 TeV LHC to completely cover all natural regions of the spectrum. Notice 
that m^o > mz/2 and m ai > rrih/2 in the spectrum which follows from constraints on the invisible 
decay width of the Z-boson and the Higgs, respectively. Finally it is interesting to focus on the 
least-tuned region which is colored in green in the spectrum. As expected for small total tuning, 
the colored sector tends to be generically lighter though there also light colored sparticles in the 
more tuned regions (colored in red and blue). 
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Figure 7: Scatter plots of the Higgs signal strengths i? 77 versus (a) Rzz and (b) R b i. Here Higgs 
production is exclusively through gluon fusion except for the bb final state where production is in 
association with a vector boson. The green, blue and red points correspond to a combined tuning 
better than 5%, between 1% and 5%, and worse than 1%, respectively. The purple points correspond 
to a gravitino LSP and have a total tuning better than 1%. The scatter plots of the couplings C 77 
versus (c) Czz and (d) C d g contain only points from the region with enhanced diphoton decay, 
1.40 < i? 77 < 1.99 (cf. Eq. (5.3)). Here the VEV fine-tuning is better than 5%. All points satisfy 
the constraints discussed in sec. 5. 

6.2 Higgs searches at the LHC 

In light of the recently discovered bosonic resonance at the LHC, the particle spectrum that we 
find leads to implications for the Higgs couplings and decays. In particular, we will consider decays 
into diphotons, dibosons (WW* and ZZ*), bb in association with a vector boson and r-f via vector 
boson fusion, which are the most relevant channels for the Higgs mass range rrih £ [124, 127] GeV. 
We write the normalized Higgs eigenstate h in terms of the gauge components h u , hj and s as 
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h = Vh,h u h u + Vh,h d hd + Vh yS s- Unless stated otherwise, we consider Higgs production through 
gluon fusion, which tends to dominate in the case of a sizable h u admixture in h. The signal 
strength to diphotons versus the signal strength to dibosons is displayed in Fig. 7a, and versus the 
signal strength to bb in association with a vector boson in Fig. 7b. First of all, let us stress that 
we can always accommodate a signal strength consistent with a SM Higgs in all channels with a 
total tuning better than 1%. However, it is also possible to explain the anomalous (at the 2cr-level) 
signal strength in the diphoton channel measured at both ATLAS and CMS. As seen in the figures, 
it is clear that an enhancement in the diphoton signal strength, consistent with the range given 
in Eq. (5.3), requires at least some reduction in the b-quark signal strength R bb . This reduction 
results from a smaller hd component in h with respect to hsM and thus also leads to a reduction 
of R rf which is, however, experimentally not very constrained at the moment. More importantly, 
as can be seen from Figs. 7a and 7b, the reduction in R bb gives rise to an increment in the diboson 
signal so that Rzz ~ R-y-y- This could be ruled out in the near future. 

However, in the regions of large total-tuning (red points which also correspond to large values 
of A) one can obtain an enhancement due to charginos [64] in the diphoton channel while keeping 
Rzz i$ R-y-y as shown in Fig. 7d. The enhancement is a consequence of large values of A that lead 
to a large negative singlet component in h, Vh, s i$ —0.3 and a sizable Higgs-chargino coupling. The 
lightest chargino 10 mainly Higgsino in composition, can then lead to enhancements in the diphoton 
channel at loop level that interfere constructively with the TV-boson contribution 



g hWW A\(Tw)+A^g h ~+~ + 



A h 1 (rw) + lytA h 1/2 (r t ) + 



(6.1) 



/2 

where t% = m^/(4mf), guww is the scale- invariant NMSSM Higgs coupling to TV-bosons and 
9h x +y- ~ V2\V 8j h is the coupling to the lightest chargino pair which being negative and large due 
to V s h and A, interferes constructively (though it is a fermion loop) with the TV-boson contribution. 

/2 



The functions A\{t) and Ay 2 (r) from loop integration are given by 



Af (r) = -2[2r 2 + 3r + 3(2r - l)/(r)]r- 2 and A h 1/2 (r) = 2[r + (r - 1)/(t)]t" 2 , (6.2) 

where /(r) = ArcSinf-y/r] 2 for r < 1, and we have used that Ay 2 (r^±) « 2. Notice that for 
V St h to be significant, one needs large values of A and thus one expects that this type of diphoton 
enhancement can only be accomplished in the large fine-tuned regions of parameter space. This is 
what we mostly found and is depicted by the red points with 1.4 < i? 77 < 1.5 shown in Fig. 7a. 

However, there are a few additional points with i? 77 > 1.7, 1.1 < Rzz % 1-32, 0.35 < R bb < 0.6 
and with total tuning better than 1% as can be seen in Figs. 7a and 7b. Most of these points are 
characterized by having large values of tan f3 > 4 and intermediate values of A ~ 1.2 such that their 
corresponding tree-level Higgs masses are close to 126 GeV and therefore their combined tuning is 
reduced as discussed at the end of sec. 5. More importantly, the lightest charginos associated with 
this region are in the range m-± G [110, 120] GeV and V^ s < —0.5 due to a reduction on Vh,h d f° r 
large tan (3, providing a large loop contribution to the diphoton decay. It is interesting that as tan j3 
becomes large, the hd component of the Higgs diminishes in magnitude, allowing for a large singlet 



°The lightest chargino is not required to be very light due to the large coupling A. 
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component and moreover, due to the T-parameter constraints, we are led to regions of moderate 
values of A where the total tuning is reduced. So what we learn from the enhanced diphoton region 
is that a sizable singlet component in the Higgs can be obtained either by: 

1. Large values of A and thus a tuned region (upper-right red corner with i? 77 ~ 1.4 in Fig. 7a). 
However this also leads to slightly heavy charginos m~± ~ fi making the diphoton enhance- 
ment mild. 

2. Large values of tan /3 which due to T-parameter constraints imply also intermediate values of 
A and thus smaller tuning (points with i? 77 > 1.6 in Fig. 7a). Since A is not as large one also 
expects that charginos can be relatively light thus providing a larger loop contribution to the 
diphoton decay. 

Furthermore most of these points correspond to a spectrum where x? is either purely singlino or 
Higgsino with m^o £ [70, 120] GeV, the xf is purely Higgsino with m~± € [110, 120] GeV and 

X2 is either purely Higgsino or singlino with m^o £ [180, 450] GeV. Thus the neutralino-chargino 
constraints discussed in sec. 5 do not apply since they are based on the assumption of a mostly 
electroweak-gaugino neutralino-chargino sector that is light. These results suggest that there is 
region of parameter space with a small amount of tuning and enhanced diphoton decay which 
deserves further study. 

The large singlet component also modifies the coupling to Z and VF-bosons and down-type 
quarks, but in the region of enhancement the down-type coupling is further suppressed than the 
W ,Z-bosoia coupling, making the production almost close to SM-like. 11 This raises all signal 
strengths due to the small BR(h — > bb) but the large singlet component allows a milder enhancement 
in the ZZ coupling compared to the diphoton coupling (which is more enhanced due to the chargino 
loop). This is shown in Fig. 7c where the ratios of Z-boson and diphoton couplings to the SM 
couplings are plotted with the result that Czz ^ C 77 , where 

Cxx = 9hxx 1 9h SM xx ■ (6.3) 

Notice that the points corresponding to this region possess a negative value of C dd ~ as shown in 
Fig. 7d, corresponding to a small but negative h d component in h, V/ l) /j d < 0. This provides an 
enhancement in the diphoton channel with down-type fermions running in the loop. However due 
to the smallness of the bottom-quark Yukawa coupling such an enhancement is insignificant and 
therefore we believe that the fact that these points correspond to a negative hd Higgs component 
is only an artifact of the scanning. 

6.3 SUSY searches at the LHC 

The naturalness constraints in the scale-invariant NMSSM have definite implications for SUSY 
searches at the LHC. We use the combined tuning in the Higgs VEV and the Higgs mass to place 
naturalness bounds on sparticle masses. We focus on the region with combined tuning better than 
1% (points colored in blue and green) and highlight the optimal region in our scan (corresponding 
to the green points) with a combined tuning better than 5%. In addition, all points have a Higgs 

11 The increase in the singlet component of h comes with a decrease in hd, while the h u component is not much 
suppressed. 
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Figure 8: Scatter plots of (a) the combined tuning T^Ti" and (b) the lightest stop mass m ( - , as a 
function of the gluino mass nig. The scatter plots in (c) and (d) show the lightest neutralino mass 
m^,o as a function of (c) the gluino mass m g and (d) the lightest stop mass m^. The green, blue 
and red points correspond to a combined tuning (T> h T, v ) better than 5%, between 1% and 5%, and 
worse than 1%, respectively. All points satisfy the constraints discussed in sec. 5. 



VEV tuning that is better than 5%. As can be seen from the spectrum shown in Fig. 6, we find 
points in our scan with a lightest stop mass up to 1.2(2.6) TeV for a total tuning better than 
5%(1%). The heavier stop can have masses up to 1.8(3.5) TeV, whereas the sbottom masses can 
reach up to 1.5(4.3) TeV for the light sbottom and 7.5(8) TeV for the heavier sbottom. The fact 
that we do not have points with sbottom masses above ~ 8 TeV is due to the chosen range of the 
input parameters in the scan (cf. table 1). 

The dependence of the combined tuning on the gluino mass is shown in Fig. 8a. As expected, 
it grows with larger gluino masses. For a combined tuning better than 5%(1%), we find gluino 
masses up to 3.0(4.6) TeV. The distribution of the lightest stop mass m£ versus the gluino mass 
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m,g is plotted in Fig. 8b. The lower bound on mi, which grows with rrig, can be understood from 
the requirement in Eq. (3.1) that there is no additional tuning in the squark mass parameters 
(cf. footnote 4). As expected, the parameter region with small stop and gluino masses is preferred 
by small fine-tuning. In Fig. 8c and Fig. 8d, we show the results of our scan in the plane m^o 
versus rrig and versus m; which is how the experimental results for gluino and stop searches are 
typically presented. The lower left corner in the plots in Fig. 8 is already excluded by stop and 
gluino searches at the LHC, which lead to the limits given in Eqs. (5.5) and (5.7). However, as stops 
and gluinos can be quite heavy for a combined tuning better than 1%, we will have to wait for the 
14 TeV LHC to completely probe the stop and gluino masses in this natural region of parameter 
space. 

The lightest electroweak charginos and neutralinos, on the other hand, are necessarily light 
because their mass is directly tied to the parameters entering the Higgs mass. In our scan, we find 
masses for the lightest chargino up to 405(625) GeV for a combined tuning better than 5%(1%), 
as well as neutralino masses up to 400(605) GeV for the lightest neutralino and 420(675) GeV 
for the next-to-lightest neutralino. Exploring this mass region for the electroweak charginos and 
neutralinos will constrain naturalness of the NMSSM at the percent level. 

Let us finally mention that for the case A mess = 100 TeV (A mess = 1000 TeV) with combined 
tuning better than 1%, gluinos are restricted to lie below 3(2.2) TeV, the lightest stop must lie 
below 2.5(1.3) TeV, whereas a light chargino m~± > 110 GeV remains possible in the spectrum. 

7 Conclusion 

A well-motivated way to accommodate the recent LHC discovery of a Higgs-like resonance at 
126 GeV is to extend the minimal supersymmetric standard model with a singlet superfield S. 
In this NMSSM, the Higgs-singlet coupling A can increase the Higgs mass already at tree-level. 
In addition, the sensitivity of the electroweak scale with respect to changes in the soft masses is 
suppressed at large A, further relieving the fine-tuning. This seems to allow the stop and gluino 
masses to be increased for the same level of tuning as in the MSSM. However, increasing the 
coupling above A > 1 typically also causes the Higgs mass to become heavier than 126 GeV. This 
then requires a further tuning among the various contributions to the Higgs mass to keep the Higgs 
light. This additional tuning limits raising the stop and gluino masses. On the other hand as A 
becomes smaller and we approach the MSSM limit, the tuning increases because we need to rely 
on large stop soft-mass parameters to raise the Higgs mass. Thus the tuning is minimized when 
A ~ 1. 

Identifying the new resonance with the lightest CP-even state in the (scale-invariant) NMSSM, 
then leads to interesting constraints on the sparticle spectrum. The naturalness constraint, which 
clearly requires light superpartners, has to confront the experimental searches from the LHC and 
electroweak precision tests. In a optimal-case scenario, with split families (including decoupled 
third-generation sleptons) and a low messenger scale of 20 TeV, we have performed a numerical 
scan over the parameter space of the (scale-invariant) NMSSM. For a tuning that is better than 
5%(1%), we find the lightest stop mass to be below 1.2(2.6) TeV, the gluino mass to be below 3.0(4.6) 
TeV and the lightest charginos and neutralinos to be below 400(600) GeV. If the messenger scale is 
increased to 100(1000) TeV, then the total tuning worsens and becomes below the 5%(2-3%) level. 

The measured DM relic abundance can be accounted for by a neutralino LSP, which is dom- 
inantly bino-like. The dominant annihilation channels are then via the Higgs and the lightest 
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pseudoscalar. In our scan, however, we mostly find Higgsino- and singlino-like neutralino LSPs, 
because the parameters determining their mass are related to the recent 126 GeV Higgs mass 
measurement. 

We also find a few points that can explain the experimentally observed enhanced diphoton 
decay rate of the Higgs. In particular, this can be achieved via a large value of A and consequently 
a sizable singlet admixture, leading to a suppression of the coupling hbb and an enhancement of 
the chargino contribution to the h — > 77 amplitude. The chargino can be relatively heavy in this 
region, since its coupling to the Higgs is enhanced by A. The more interesting possibility occurs 
for larger values of tan f3, implying intermediate values of A and small fine-tuning. In this region, 
the lightest chargino has to be light to give a substantial enhancement of h — > 77. Indeed, we find 
chargino masses of 110 — 120 GeV, which will be probed experimentally soon. We believe that this 
region of small fine-tuning and enhanced Higgs diphoton rate deserves further study. 

In summary, our analysis shows that in the optimal-case scenario of a low messenger scale (20 
TeV) and light colored third-generation sparticles there are still sizable regions in the parameter 
space of the scale-invariant NMSSM where the tuning is better than 5%. This still allows for stop, 
gluino, chargino and neutralino masses above the current LHC limits that will be eagerly searched 
for at the 14 TeV LHC, thereby probing naturalness down to the percent level. 
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Appendix 

A Mass matrices 

In this appendix, we collect the mass matrices of stops, sbottoms and the Higgs-singlet sector 
(except for the CP-even states for which the mass matrix has been discussed already in sec. 2). 

In the CP-odd sector of the Higgs-singlet scalars, we first decouple the Nambu-Goldstone boson 
"eaten" by the massive Z-boson. The remaining CP-odd state in the Higgs sector is given by the 
linear combination A = sin f3 h^j + cos f3h u j. The mass matrix in the basis Acp- Q dd = (A,si) 
then takes the form 



Diagonalizing the mass matrix, we obtain the mass eigenstates a± and (12, where we label the states 
such that their masses satisfy m ai < m a2 . 



CP-odd 




(A.l) 
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Similar to the MSSM, gauginos and Higgsinos mix after electroweak symmetry breaking. In 
addition, there is a new mixing mass term proportional to A between the singlino and Higgsinos. 
In the gauge-eigenbasis ip° = (B, W 3 , H®, H®, S), the neutralino mass matrix is then given by 



/ Mi — cos j3 sin 9w m z sin (3 sin 6wmz 
M2 cos (3 cos Owmz — sin f3 cos Owmz 

—fj, 



v . . 



\ 


-Xv sin /3 
-Xv cos /3 



(A.2) 



The charginos mix as in the MSSM. In the gauge-eigenbasis V ;=t 
mass matrix is therefore given by the usual 2x2 block form 



(W+,H+,W-,H7), the 



M, 



X T 
X 



where X 



M 2 

V2mw cos /3 



\^2mw sin (3 
/i 



(A.3) 



and the superscript T denotes the transpose. 

The stops and sbottoms also mix as in the MSSM. The mass matrix of the stops is given by 



f a t 


bt\ 




V bt 


ct J 


' ( 



™t + (2 " I 



sin 2 9w) fn 2 z cos 2 P 



m t (A t 



tan /? 



m 



us 



+ mf + I sin 2 9wm 2 z cos 2/3 / ' 

(A.4) 



where the ^-term is defined such that AC so jt = ytAtu^Q^H u +c.c.. Diagonalizing the mass matrix, 
we obtain the mass eigenstates t\ and t 2 with masses given by 



m l,ti = l(at + c t ±A t ), 



(A.5) 



where At = \J (at — ct) 2 + 46 2 . Similarly, the mass matrix of the sbottoms is given by 



a b b b 
h c b 



rriQ 3 + (~\ + \ sin 2 W ) m| cos 2/3 



rrib(A b — /xtan/3) 
m d 3 + m b ~ I sm2 Owm 2 ? cos 2/3 



(A.6) 



and the masses of the mass eigenstates b\ and b 2 follow from Eq. (A.5) with the replacement 
(a t ,bt,c t ) ->■ (a b ,b b ,c b ). 



B Effective Potential 

In this appendix, we give the dominant contributions to the Coleman- Weinberg potential, due to 
third- generation (s) quarks and the Higgs-singlet sector itself. Defining the function 

f( m )^m 4 (ln^-l) , (B.l) 
V ™ so ft 2 J 

we can express these contributions in a compact form. 
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The contribution due to third-generation (s)quarks to the effective potential is given by 



Vi. 



64vr 2 



2 E f( m + 2 E K) - 4 ^ m *) - 4 ^ m <>) 



i=l 



i=l 



(B.2) 



where iV c = 3 and mj. denotes the stop masses (given in Eq. (A. 5)) and the sbottom masses. 
The top and bottom mass are given by m t = ytv u and m& = ybVd, respectively. All masses are 
evaluated at the renormalization scale m so ft = ^/mQ 3 m U3 . The different prefactors (2 for sfermions 
and -4 for fermions) originate from the degrees of freedom and the different statistics. 
The contribution due to the Higgs-singlet sector (including superpartners) is given by 



Vi 



1 



l,Higgs 



64tH 



3 

£ 

i=i 



f(m s 



2 

+E 

i=l 



2 E/( 



«=i 



4E/K f ; 



8=1 



(B.3) 



The masses are functions of the parameters in the Higgs sector, especially of A and k, which can 
become sizable. An explicit calculation of this contribution to the effective potential is tedious, 
since it involves the diagonalization of a 3 x 3 and a 5 x 5 mass matrix. 

The calculation of the minimization conditions and the fine-tuning measure requires the knowl- 
edge of several derivatives of the mass eigenvalues with respect to the VEVs as well as parameters 
in the Lagrangian. The existence of an analytic expression for the masses is not necessary, however, 
since an expression for the derivatives can be obtained by differentiating the characteristic polyno- 
mial and solving the resulting linear equation for the derivative. This is similar to the calculation 
of the fine-tuning measure discussed in sec. 4. 



References 

[1] S. P. Martin (1997), hep-ph/9709356. 

[2] U. Ellwanger, C. Hugonie, and A. M. Teixeira, Phys.Rept. 496, 1 (2010), 0910.1785. 

[3] R. Barbieri, L. J. Hall, Y. Nomura, and V. S. Rychkov, Phys.Rev. D75, 035007 (2007), 
hep-ph/0607332. 

[4] L. J. Hall, D. Pinner, and J. T. Ruderman, JHEP 1204, 131 (2012), 1112.2703. 
[5] S. Abel, Nucl.Phys. B480, 55 (1996), hep-ph/9609323. 

[6] C. Panagiotakopoulos and K. Tamvakis, Phys.Lett. B446, 224 (1999), hep-ph/9809475. 
[7] T. Gherghetta and A. Pomarol, Phys.Rev. D67, 085018 (2003), hep-ph/0302001. 
[8] R. Sundrum, JHEP 1101, 062 (2011), 0909.5430. 

[9] T. Gherghetta, B. von Harling, and N. Setzer, JHEP 1107, 011 (2011), 1104.3171. 

[10] R. Harnik, G. D. Kribs, D. T. Larson, and H. Murayama, Phys.Rev. D70, 015002 (2004), 
hep-ph/0311349. 

[11] E. Hardy, J. March-Russell, and J. Unwin, JHEP 1210, 072 (2012), 1207.1435. 



32 



[12] R. Franceschini and S. Gori, JHEP 1105, 084 (2011), 1005.1070. 
[13] Z. Kang, J. Li, and T. Li, JHEP 1211, 024 (2012), 1201.5305. 

[14] S. King, M. Muhlleitner, and R. Nevzorov, Nucl.Phys. B860, 207 (2012), 1201.2671; 
S. F. King, M. Muhlleitner, R. Nevzorov, and K. Walz (2012), 1211.5074. 

[15] D. A. Vasquez, G. Belanger, C. Boehm, J. Da Silva, P. Richardson, et al., Phys.Rev. D86, 
035023 (2012), 1203.3446. 

[16] J.-J. Cao, Z.-X. Heng, J. M. Yang, Y.-M. Zhang, and J.-Y. Zhu, JHEP 1203, 086 (2012), 
1202.5821; 

J. Cao, Z. Heng, J. M. Yang, and J. Zhu, JHEP 1210, 079 (2012), 1207.3698. 

[17] K. Agashe, Y. Cui, and R. Franceschini (2012), 1209.2115. 

[18] S. R. Coleman and E. J. Weinberg, Phys. Rev. D7, 1888 (1973). 

[19] J. D. Hunter, Computing In Science & Engineering 9, 90 (2007). 

[20] A. Arvanitaki, N. Craig, S. Dimopoulos, and G. Villadoro (2012), 1210.0555. 

[21] U. Ellwanger, J. F. Gunion, and C. Hugonie, JHEP 0502, 066 (2005), hep-ph/0406215. 

[22] U. Ellwanger and C. Hugonie, Comput.Phys.Commun. 175, 290 (2006), hep-ph/0508022. 

[23] G. Belanger, F. Boudjema, C. Hugonie, A. Pukhov, and A. Semenov, JCAP 0509, 001 (2005), 
hep-ph/0505142. 

[24] A. Djouadi, M. Drees, U. Ellwanger, R. Godbole, C. Hugonie, et al., JHEP 0807, 002 (2008), 
0801.4321. 

[25] G. Belanger, F. Boudjema, P. Brun, A. Pukhov, S. Rosier-Lees, et al., Comput.Phys.Commun. 
182, 842 (2011), 1004.1092. 

[26] E. Aprile et al. (XENON100 Collaboration), Phys.Rev.Lett. 109, 181301 (2012), 1207.5988. 

[27] M. E. Peskin and T. Takeuchi, Phys.Rev. D46, 381 (1992). 

[28] R. Barbieri and G. Giudice, Nucl.Phys. B306, 63 (1988). 

[29] M. Perelstein and B. Shakya (2012), 1208.0833. 

[30] S. Cassel, D. Ghilencea, and G. Ross, Nucl.Phys. B835, 110 (2010), 1001.3884. 

[31] G. Larsen, Y. Nomura, and H. L. Roberts, JHEP 1206, 032 (2012), 1202.6339. 

[32] S. Heinemeyer, W. Hollik, and G. Weiglein, Comput.Phys.Commun. 124, 76 (2000), hep-ph/ 
9812320. 

[33] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B (2012), 1207.7214. 

[34] S. Chatrchyan et al. (CMS Collaboration), Phys.Lett. B716, 30 (2012), 1207.7235. 

33 



[35] Tech. Rep. ATLAS-CONF-2012-162, CERN, Geneva (2012). 

[36] Tech. Rep. ATLAS-CONF-2012-170, CERN, Geneva (2012). 

[37] Tech. Rep. CMS-HIG- 12-045, CERN, Geneva (2012). 

[38] Tech. Rep. CMS-PAS-HIG-1 1-024, CERN, Geneva (2011). 

[39] Tech. Rep. CMS-PAS-HIG-12-041, CERN, Geneva (2012). 

[40] J. Beringer et al. (Particle Data Group), Phys.Rev. D86, 010001 (2012), URL http://link. 
aps . org/doi/10 . 1103/PhysRevD . 86 . 010001. 

[41] Tech. Rep. ATLAS-CONF-2012-072, CERN, Geneva (2012). 

[42] G. Aad et al. (ATLAS Collaboration), Phys.Lett. B715, 44 (2012), 1204.6736. 

[43] J. Barnard, B. Farmer, T. Gherghetta, and M. White (2012), 1208.6062. 

[44] Tech. Rep. ATLAS-CONF-2012-145, CERN, Geneva (2012). 

[45] Tech. Rep. ATLAS-CONF-2012-151, CERN, Geneva (2012). 

[46] Tech. Rep. SUS-12-028-pas, CERN, Geneva (2012). 

[47] Tech. Rep. ATLAS-CONF-2012-106, CERN, Geneva (2012). 

[48] G. Aad et al. (ATLAS Collaboration) (2012), 1208.1447. 

[49] G. Aad et al. (ATLAS Collaboration) (2012), 1208.2590. 

[50] Tech. Rep. CMS-PAS-SUS-12-023, CERN, Geneva (2012). 

[51] Tech. Rep. ATLAS-CONF-2012-154, CERN, Geneva (2012). 

[52] S. Chatrchyan et al. (CMS Collaboration), JHEP (2012), 1209.6620. 

[53] Tech. Rep. CMS-PAS-SUS-12-022, CERN, Geneva (2012). 

[54] G.-C. Cho and K. Hagiwara, Nucl.Phys. B574, 623 (2000), hep-ph/9912260. 

[55] K. Hagiwara, S. Matsumoto, D. Haidt, and C. Kim, Z.Phys. C64, 559 (1994), hep-ph/ 
9409380. 

[56] S. P. Martin, K. Tobe, and J. D. Wells, Phys.Rev. D71, 073014 (2005), hep-ph/0412424. 

[57] Y. Amhis et al. (Heavy Flavor Averaging Group) (2012), 1207.1158. 

[58] LHCb collaboration (2012), 1211.2674. 

[59] F. Domingo and U. Ellwanger, JHEP 0712, 090 (2007), 0710.3714. 

[60] D. Larson, J. Dunkley, G. Hinshaw, E. Komatsu, M. Nolta, et al., Astrophys.J.Suppl. 192, 16 
(2011), 1001.4635. 

34 



[61] A. de Gouvea, T. Moroi, and H. Murayama, Phys.Rev. D56, 1281 (1997), hep-ph/9701244. 

[62] G. Giudice and R. Rattazzi, Phys.Rept. 322, 419 (1999), hep-ph/9801271. 

[63] N. Fornengo, A. Riotto, and S. Scopel, Phys.Rev. D67, 023514 (2003), hep-ph/0208072. 

[64] K. Schmidt-Hoberg and F. Staub, JHEP 1210, 195 (2012), 1208.1683. 



35 



